Low Temperature Upper Critical Field Anomalies in Clean Superconductors 
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We interpret the upper-critical-field anomalies observed in some high- temperature superconduc- 
tors as resulting from the proximity to a zero-temperature quantum critical point. We estimate the 
shape of the phase boundary between the normal and the superconducting phase by modeling the 
zero temperature critical point as the second-order endpoint of the first-order melting line of the 
vortex lattice. 
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In two recent letters, Mackenzie et al. and Osofsky et. al. observed that in overdoped high temperature 
superconductors the upper critical field can be measured at very low temperature and displays a very steep rise as 
the temperature is decreased Mackenzie et. al. studied the single layer system Tl2CuOe+s which is overdoped by 
incorporating excess oxygen, while Osofsky et al. studied Bi2SrCu02- 

The positive curvature of the Hc2 vs. T curve and its rapid incr anomalous low temperature cannot be explained 
by the classical WHH (Werthamer Helfand Hohenberg) theory. In fact, these anomalies have triggered a large number 
of theoretical interpretations. Schofield and Wheatley suggested that these anomalies are a manifestation result from 
the Luttinger liquid behavior in the normal state . Alexandrov, Bratkovsky and Mott |^ account for the anomalous 
curvature of the Hc2VsT curves in terms of bipolaron superconductivity and Brandow suggested that the negative 
curvature is the result of pair breaking Mackenzie et. al. |^ have suggested that the upper critical field anomalies 
are related to a strong temperature dependence of the effective mass. 

In this paper we suggest a very different origin of the upper critical field: the proximity to a zero temperature critical 
point. We first argue, on very general grounds, that the observed low-temperature behavior implies the existence of 
a thermodynamic singularity at zero temperature which we characterize in terms of two critical exponents. Then we 
propose a simple model for the relevant combination of exponents controlling the shape of the low temperature upper 
critical field line which agrees well with the experimental observations. 

The microscopic considerations involve the melting of the vortex lattice, a problem which has been a subject of 
intense study Most of the work in this area has concentrated on the classical statistical mechanical aspects of 
this problem with the exception of the recent work by Blatter and collaborators which carried out a microscopic 
calculation of the quantum and thermal fluctuation of the vortex lattice to one loop order. 

Using standard thermodynamic identities, one can relate the slope of the upper critical field curve separating the 
normal and the superconducting phase to the change in entropy and magnetization across the transition line: 

Sn ~ Ss _ dHc2 N 

M„ - M. " ~dT~ ^ ' 

If the superconductor-to-metal transition is of the first order, the latent heat and the magnetization jump are finite. 
If it is second order, the left hand side of equation (|l|) should be understood as a derivative. At zero temperature 
^(r ~ 0, M) = On the normal and the superconducting phase. If S is a regular function as M approaches the upper 
critical field line and the temperature tends to zero, then Eq. |l| implies that the slope of the upper critical field curve 
vanishes at zero temperature. In the experiments of Ref. 0] and the slope of the upper critical field curves diverges 
rather than vanishes as the temperature approaches zero, implying the existence of a singularity in the free energy at 
zero temperature, namely a quantum critical point. 

Below the upper critical dimension, a second-order transition from an Abrikosov type // superconductor to a normal 
metal as a function of field is parametrized by two independent exponents v and z which control the divergence of a 
length scale and of a time scale as the critical point is approached. The free energy per unit volume has a singular 
part above and below the superconducting transition which behaves as 

f = A{H- H,2r^''+'^ (2) 

At finite temperature scaling implies 

/(T, H) = A{H - H,2{T = Oyy ^^+''^g{T[H - H,2{T - 0)]^^) (3) 

The location of the upper critical field vs. temperature line is given by the position in the temperature-field plane 
where the free energy is singular. If we denote the singular point of the scaling function g by Xc, we obtain a connection 
between the shape of the upper critical field at low temperatures and the critical exponents: 

Hc2 (T) = H,2{T = 0) ~ XcT^ . (4) 

The steep decrease in Hc2 observed experimentally is then a measure of the product of the dynamical critical exponent 
and the correlation length exponent of the zero- temperature critical point. The scaling assumption connects the shape 
of the critical line to singularities in other physical quantities. We expect the linear term in the the normal-state 
specific heat to have a singular part behaving as 7 w {H — Hc2y^'^~^^ while the a.c. susceptibility acquires a singularity 
of the form x ~ {H — Hc2Y^-'''~^^^~^ as we approach the upper critical field at low temperatures. 

To estimate the exponent which determines the shape of the upper critical line at low temperatures we need a more 
detailed picture of the critical point. Presently there is no microscopic theory of the superconductor to metal phase 
transition in the presence of a magnetic field in three dimensions. To make progress we assume that at any finite 
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temperature, the superconductor to normal metal transition is ( weakly) first order, an assumption which is supported 
by some renormalization group calculations [Tl| ] and that this first order line ends at zero temperature in a quantum 
critical point, whose existence is strongly suggested by the experimental data. Then at any finite temperature, we 
are dealing with the first order melting transition of the vortex lattice, and the proximity to the second order zero 
temperature critical point is taken into account by using a renormalized parameters in the determination of the 
melting line. 

More precisely we regard the phase transition at finite temperatures as a (weakly first order) melting of the Abrikosov 
lattice of an anisotropic three-dimensional superconductor. Then we estimate the locus of the transition line by a 
modified version of the Lindemann criterion [|l5| that takes into account the zero-temperature critical behavior near 
iJc2(0). Our basic idea is that at any finite temperature we can integrate out the quantum fluctuations to obtain an 
effective action that describes the finite temperature transition. After this integration of the quantum fluctuations 
is carried out the proximity to the quantum critical point is contained solely in two (renormalized by quantum 
fluctuations ) parameters: a renormalized stiffness and a renormalized Lindemann number. 

We assume that melting occurs when the root mean square displacement of a vortex v<u^> due to thermal 
fluctuations becomes of the order of dy a quantity related to the distance between the normal regions surrounding 
nearby vortices. The square of the distance between the centers of the vortices 1{H) deflned by {^^j^ = ^), is set 

by the external magnetic field. The area covered by normal region surrounding the core is given by with 
Ic = l{Hc2{T)) . The quantity = P — vanishes as H approaches Hc2'- 

dy~l ^^^(^^ . (5) 



Our modification of the Lindeman criteria is based on a the following picture. After time averaging over many 
osillations of the vortex cores, each vortex is associated with an area ^^^^ — h const < > of normal phase. Melting 
occurs when these normal areas cover the whole sample, that is: 



V< > = ady (6) 

Here a is a dimensionless constant similar in spirit to the Lindemann number. Based on this analogy we expect it to 
be smaller than one. 

We stress that this is our main assumption and is of a phenomenological nature. In this equation Hc2{T) is the upper 
critical field computed without taking into account the thermal fluctuations which are responsible for the melting of 
the vortex lattice. 

The thermal contribution to the root-mean-square displacement \/< >, is proportional to the temperature and 
inversely proportional to a typical elastic constant. For the Abrikosov lattice this estimate is complicated by the fact 
that some elastic moduli are highly non local and softer at short wavelengths than at long wavelengths. At zero 
wave vector the bulk modulus cii(O) and the tilt modulus 044(0) are non critical while the elastic shear modulus is 
given by C66(0) ~ (i?c2 — H)^. The expected softening of the bulk and tilt moduli does occur at large wavevectors.For 
kh with kh ~ {Hc2 — H)^ ^ C44{q) « (fc^/g)^C44(0), cii(q) « (fc/i/(?)^cii(0) giving rise to strong infrared divergences 
in the evaluation of •\/< >. Fortunately the relevant estimate of the vortex lattice mean displacement taking into 
account the nonlocality of the bulk moduli has been carried out by Brandt and by Houghton et. al. in the 
anisotropic case. They showed that as H approaches Hc2 the most divergent contribution to the thermal displacement 
has the form < B,P n\T^^^)ikT with B, ^ (-^^j^^)^(^)^.^. 

Here M and Mz are the masses in the Landau-Ginzburg Hamiltonian, c(0) denote the elastic constants at q = 
, A is an ultraviolet cutoff of the order of the vortex lattice zone boundary wave vector, and k is the ratio of the 
penetration depth to the coherence length. Inserting the field dependence of the g = elastic moduli they obtained 

<u^>^Bl^>e(^\ (7) 

where B is a constant equal to 2.26 x \Q~^K~'^ . 

To obtain the melting curve, which we denote H„i{T), we insert Eqs. (^ and ( |^) into the generalized Lindemann 
equation (||). 

The main difference between our equation and the usual Lindemann criterion which takes dy to be a non-critical 
constant of the order of the magnetic length, is the critical field dependence of dy and < >. 
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Our assumption for (||) which in a more general case could involve an arbitrary exponent is akin to a mean-field 
approximation incorporating the physics of the proximity to the zero-temperature critical point. Temperature appears 
explicitly in equations ^and|^ and also implicitly in Hc2{T). This temperature dependence only introduces analytic 
(proportional to (T/Tc)'^) corrections which are negligible compared with the non-analytic terms that we derive below. 
We therefore evaluate Hc2{T) at zero temperature and obtain an equation for Hm{T) given the value of Hc2{T — 0). 

A plot of (5) together with the data of Mackenzie et. al. in a range including 2 decades of reduced temperature 
is shown in Fig. 1. The theoretical fit used the values Hc2{0) = 17 AT and T* = 3.99K We find that this estimate 
fits the low-temperature portion of the data of Ref. ( 0] ) remarkably well. At higher temperatures and lower fields 
the curves depart from this simple power law behavior, but we do not expect our considerations to apply far from the 
zero-temperature critical point. At the lowest temperatures measured, there is a hint of a crossover to a power law with 
a power much closer to one. More experimental data with smaller error bars in the low temperature region is needed 
to determine wether our theory holds in the very low temperature region or wether there is a relevant perturbation 
such as disorder which is responsible for another type of critical behavior.. The theory has a free parameter which 
is the value of a in Eq. (|). This is very similar to the parameter c in the Lindemann theory of melting [Tsl . It 
determines the characteristic temperature T* via T* « 't%ii"°i/" ■ Using the experimental values for the anisotropy 

ratio in Tl2Ba2CuOQ+s, ~ 31.6 and a Ginzburg parameter n w 200 iQ, we find that our value of T* 
corresponds to a ~ .1. 

The zero point motion gives of course a T-independent contribution to the mean-square displacement, which can be 
added to the right hand side of the Lindemann criterion. If it is only weakly field dependent it can be considered as 
merely a renormalization of the core radius, or effectively a renormalization of Hc2 at T = 0. The one-loop calculation 
of Blatter et al. shows that this is the case We stress that once quantum corrections are explicitly included in a 
renormalized value of i?c2(0), the elastic moduli have to be recalculated in a consistent fashion so that they vanish 
when the field equals the renormalized value of -ffc2- 

Other scenarios, such as the naive application of the ordinary Lindemann criterion, or the theory of Kosterlitz- 
Thouless two- dimensional melting give an exponent of 2/3 and do not agree well with the experimental data. 

We now turn to the reasons why the upper-critical-field singularity was observed only in Refs. |^] and |^ and not 
in other high temperature superconductors where fluctuations are clearly important. First we note the the size of the 
upper critical field is not anomalous if one takes into account the short coherence length of these systems. In fact 
using the back- of-the-envelope estimate for the zero temperature upper critical field (po = Bn^^ and an estimate of 
^ w 60A" we obtain B « 16T. For our considerations to be applicable a sharp phase transition should take place 
between the superconducting and normal phase without an intermediate wide crossover region, the so called "vortex 
liquid regime" . In the materials we discuss in this paper sharp resistive transitions take place. The sharpness is 
due to a combination of the purity of the samples and their strong anisotropy which eliminates the kinetic barriers 
responsible for the vortex liquid phase. Disorder and the smaller anisotropy of YBCO cause pinning and vortex 
entanglement which result in a new intermediate asymptotic regime, the vortex liquid, broadening the transition 
between the superconducting and the normal phase. 

Our ideas relating the Hc2 anomalies to a zero-temperature quantum critical point can be tested experimentally by 
looking for critical behavior in specific heat and susceptibility measurements. 

If the quantum critical point is connected to melting, we expect that controlled addition of impurities to the sample 
will result in a broadening of the transition. 

The considerations in this paper are largely phenomenological and motivated by experiments. To justify these 
ideas from microscopic considerations, one is led to the problem of the quantum melting of the vortex lattice. If the 
considerations presented here are correct, quantum effects can turn the weakly first order finite temperature melting 
transition into a continuous transition at zero temperature. 
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FIG. 1. Upper critical field in Tcsla vs temperature. The continuous line is given by eq. 7 with Hc2(Q) = 17.4T and 
T* = 3.99 while the stars are the experimental points from Ref. [1] 



5 



